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Global results on reset-induced periodic trajectories of planar systems 


Andrea Bisoffi^ Fulvio Forni^, Mauro Da Lio' and Luca Zaccarian^ 


Abstract —We study the existence of asymptotically stable 
periodic trajectories induced by reset feedback. The analysis 
is developed for a planar system. Casting the problem into 
the hybrid setting, we show that a periodic orbit arises from 
the balance between the energy dissipated during flows and the 
energy restored by resets, at jumps. The stability of the periodic 
orbit is studied with hybrid Lyapunov tools. The satisfaction 
of the so-called hybrid basic conditions ensures the robustness 
of the asymptotic stability. Extensions of the approach to more 
general mechanical systems are discussed. 

I. Introduction 

Starting from the important theorem of Poincare- 
Bendixson, many theoretical efforts have been made in the 
characterization of periodic orbits for planar continuous-time 
nonlinear systems, motivated by the pervasive presence of 
oscillators in electronics, mechanics and biology [5], [3]. A 
recent research direction seeks to extend this effort to the 
hybrid setting, namely to the context where, for a planar 
dynamical system, a suitable interplay of continuous flow 
and discrete jumps of the solutions leads to the existence 
of attractive periodic hybrid trajectories. The relevance of 
this topic in engineering is readily shown by the studies on 
bipedal robotic walking, where periodic hybrid trajectories 
arise from the combination of the free motion of the legs 
(continuous flow) with the impulsive action of the impacts 
at ground contact (discrete jumps) [14], [13]. 

The paper provides a stability analysis of hybrid periodic 
trajectories for planar mechanical systems based on the 
hybrid Lyapunov stability tools in [1]. The main motivation 
for the paper comes from the literature on variable impedance 
actuators, typically adopted in robotics. Strongly inspired by 
biological musculoskeletal systems, these actuators have a 
tunable stiffness and/or damping, which play a relevant role 
to improve motion efficiency [4], [9], [8], [7]. For multi¬ 
body systems with frequency separation between hrst and 
subsequent natural modes, [6] and [9] show that periodic 
oscillations can be obtained by means of simple switching 
control laws tuned only on the hrst natural mode. Taking 
advantage of a number of hybrid tools, we revisit and extend 
the results in [6]. We model the dynamics in [6] as a hybrid 
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system and we show the existence of a unique (hybrid) 
periodic orbit arising when the energy dissipated during 
flow balances the energy restored by the reset-control action 
at a jump. The stability analysis exploits hybrid Lyapunov 
methods. In particular, the asymptotic stability of the periodic 
orbit follows from the decay along system trajectories of 
a suitable Lyapunov function tailored on the kinetic and 
potential energies just after and just before a jump. 

The most relevant advantage of the approach is the in¬ 
trinsic (in-the-small) robustness of asymptotic stability [1, 
Chapter 7], which makes possible the use of the reset 
feedback law in applications. The robustness of the design 
guarantees that the stability of the attractor persists, and is 
degraded with continuity, in the presence of small parameter 
perturbation or when the instantaneous reset law is replaced 
by a (sufficiently) fast continuous actuation. 

The asymptotic stability of the attractor holds for any 
parameters conhguration that allows for a unique periodic 
orbit. This follows from the fact that the Lyapunov function 
is based on the mechanical energy just after and just before a 
jump. Its minimum is represented on the phase space by the 
set of points such that the dissipated and restored mechanical 
energy are balanced. Its decay is a natural consequence of 
the mechanical features of the system. Indeed, no explicit 
characterization of the periodic orbit is required. 

We anticipate that our Lyapunov-based analysis has simi¬ 
larities with the classical Poincare analysis of periodic orbits. 
The level sets of the Lyapunov function are univocally iden- 
tihed by the points of the hyperplane at which resets occur. 
This hyperplane plays the role of a Poincare section. Namely, 
along the portion of trajectory starting from and returning to 
this hyperplane, the overall decay of our Lyapunov function 
captures the convergence of the return map towards the 
hxed point. The advantage of a Lyapunov analysis is the 
characterization of the basin of attraction of the periodic 
orbit. In this sense, our approach is close in nature to the 
analysis of the rimless wheel in [12]. 

A promising future direction from our results is to exploit 
the hybrid framework to provide mixed continuous-discrete 
control strategies to optimize motion efficiency. 

The paper is organized as follows. The hybrid dynamics 
is discussed in Section |II] Sections [nl| and |IV] provide 
conditions for the existence of periodic hybrid trajectories 
and for their stability. Technical proofs are in the Appendix. 
Simulations in Section |V] illustrate the convergence towards 
the unique hybrid periodic orbit of the system. A comparison 
with the literature and further discussions are reported in 
Section 


II. System description 

Based on [ 6 ], consider the classical mass-spring-damper 
mechanical system 

mq-\-cq + k{q —9) =Q ( 1 ) 


with mass, damping and elastic constants respectively m, c, 
k. q is the displacement of the mass and 0 is the control input. 
The elastic force provided by the spring is proportional to 
the difference q — 9. The role of 9 is to enforce a variation 
in the stored potential energy of the spring. Following [4], 9 
could model the effect of the slow preloading of the spring 
during the flight phase of a hopping robot, which is then 
released by a clutch mechanism when touching the ground. 

In what follows 9 is piecewise constant: it switches be¬ 
tween 9 G {—9/2, 9/2} when the trajectories of the system 
pass through the hyperplane defined by {(q, g) G | q—9 — 
0 }. 0 > 0 is a design parameter corresponding to the amount 
of potential energy loaded in the spring at switches. Switches 
on 0 can be considered as the limit of a very fast continuous 
action on the spring, a kick of energy, rapidly moving 9 from 
one value to the other in {— 0 / 2 , 0 / 2 }. 

With coordinates xi := q — 9 and X 2 := q, the dynamics 
of the system can be represented according to the hybrid 
formalism in [1] as follows. Since 0 is constant, the flow 
dynamics reads 


X = f(x) 




, X G C. 


The flow set C enabling flow dynamics is given by 


( 2 a) 


C = {{xi,X2) G X 1 X 2 < 0} 

U {(xi, X 2 ) G : |xi| > 0, X 1 X 2 > 0} . (2b) 


The jump dynamics reads 


G G{x) := 


0sgn(x2) 

X2 


, X GV 


where 


sgn(x2) = 


( 2 c) 


Jsign(x2) ifx2 7^0 
[{ 1 ,- 1 } ifX 2 = 0 . 

The jump set V enabling jump dynamics is given by 
V = {(xi,X2) G xi = 0} . 


( 2 d) 



Figure [T] provides a graphical illustration of the flow and 
jump set on the system phase plane. ( | 2 d| ) guarantees that 
jumps occur when q — 9 = xi =0. For X 2 7 ^ 0, we have 
that |xi| is reset from 0 to lx/"! = 0 that is, |g+ — 0 +| = 
|g — 0+1 = 0. Indeed, the reset corresponds to a switch in 
the equilibrium position of the spring, through actuation. We 
do not reset the mass position q. 

The behavior of the solutions is illustrated in Figure 
which for a system with parameters m. = 1 kg, c = 0.3 Ns/m, 
k = 1 N/m, 0 = 0.2 m, for two different initial conditions. 
The two trajectories converge asymptotically to an attractor 
defined by the image of a hybrid periodic trajectory, where 
periodicity must be intended in a hybrid sense as clarified in 
the next section. 

Remark 1: In pcj l we used the set-valued mapping sgn to 
guarantee that the graph of the jump map x 1 — G{x) is a 
closed set. This feature ensures the outer semicontinuity of 
G. Outer semicontinuity of G combined with the continuity 
of / and with the fact that C and T) are closed sets guarantees 
that hybrid system Q satisfies the hybrid basic conditions [ 1 , 
Assumption 6.5]. They guarantee regularity of the solution 
set and robustness to small perturbations [ 2 ]. j 


Phase portrait Phase portrait 



Fig. 2: Phase plot of hybrid trajectories. Left: (xi^O; 2 : 2 , 0 ) := 
(0.1,—0.05). Right: (xi,o,X2,o) := (0.5,—0.05). 

III. Hybrid periodic orbits 

The notion of periodicity for a hybrid trajectory is a 
straightforward extension of the usual notion of periodicity. 

Definition 1: Given any hybrid system T-L := (C, F, D, G), 
a hybrid periodic trajectory is a complete solution x for 
which there exists a pair (T, J) with either T G K>o and J G 
Z>o or T G K>o and J G Z>o such that (t,J) G dom(x) 
implies (t + T, j + J) G dom(x) and, moreover, 

x(t,j) = x(t + T,j + J). (3) 

The image of x is a hybrid periodic orbit. 

The following standing assumption on the parameters of 
the hybrid dynamics Q is necessary for the existence of a 
nontriviaQ hybrid periodic trajectory. 

Assumption 1: 0, m, c and k are strictly positive. The 
roots of ms^ -f cs -f fc = 0 are complex conjugate, that is 

_ A < 0 

Assumption guarantees that mxi + cxi + kxi = 0 
is an underdamped mechanical system [10, Chap. 2.2]. 

*A nontrivial hybrid periodic orbit comprises more than one point. 































When the system is not underdamped there is no guarantee 
that a nontrivial hybrid periodic trajectory exists. With real 
eigenvalues in the flow map, the trajectories of the system 
may converge to the origin according to the direction of the 
eigenvector corresponding to the slowest eigenvalue, lying in 
the second/fourth quadrant for m,c,k > 0. In such a case, 
the solutions to Q exhibit at most one jump and the origin 
is a globally asymptotically stable equilibrium. 

The existence of a nontrivial hybrid periodic orbit follows 
from energy considerations. Consider the H-shaped curve Co 
represented in Figure given by the set 

Co = {x € : \xi I = 0, X 1 X 2 > 0} 

U {x e : |a;i| < 0, X 2 = 0}. (4) 

Under Assumption[2 the trajectories starting from Cq nec¬ 
essarily flow until they reach T). More specifically, flowing 
solutions from any a; G C in forward (respectively, backward) 
time reach set V (respectively, Co) in finite time because of 
the revolving nature of the flow trajectories. The following 
quantities are thus well defined. 

• Backward energies. Denoting by {xib, X 2 b) the intersec¬ 
tion with Co after flowing in backward time from x, 

Tb{x) ■■= ^mxlb Ub{x) := ^kxjb (5a) 

are the backward kinetic and backward potential ener¬ 
gies, respectively. 

• Forward energies. Denoting by {xif,X 2 f) the intersec¬ 
tion with T) after flowing in forward time from x, 

Tf{.x) := ^mxlf (5b) 

is the forward kinetic energy. 

Figure shows the level sets of T^, Ub and Tf, which 
correspond indeed to flowing portions of solutions to Q. 




Fig. 3; Set Co. Flowing solutions. The curves from Co to T) 
are level sets of Tf,, Ub and Tf. 

For each x f 0 the quantity Tb{x) + Ub{x) is the total 
mechanical energy of the system right after a jump. The 
quantity Tf(x) is the total mechanical energy of the system 
after a maxima]^ flow, that is, right before a jump. The 
difference between these two energies corresponds to the 
dissipation during flows. 

^In the same sense of maximal solutions in [1, Definition 2.7], namely 
that it can not be extended further. 


The reset of 6 injects energy into the system in the form 
of potential energy. This fact and the central symmetry of 
the phase portrait (namely, if </) is a solution to 
is a solution as well) imply that a hybrid periodic orbit 
corresponds to the set of points satisfying the energy balance 

Tb[x) + Ub{x) = Tf{x) + ^k9'^ (6) 

where the last term represents precisely the potential energy 
injected by a reset. Given the mentioned central symmetry, 
X f 0 belongs to a periodic hybrid orbit only if Ub{x) = 
so that (jhi is equivalent to Tb{x) = Tf{x). 

Existence and uniqueness of the hybrid periodic orbit 
follows from the monotonicity of the dissipation with respect 
to initial conditions on Cq. In fact, for x = {xi^X 2 ) G Co the 
dissipation is a strictly increasing function of |a;|. This is 
precisely stated in the next lemma, proven in the Appendix. 

Lemma 1: Consider any solution x to 0 flowing from Cq 
at ordinary time ti to V at ordinary time t 2 > fi and define 
the total mechanical energy at x as E{x) = ^mX 2 + ^kxf 
The dissipated energy E{x{ti)) — E{x{t 2 )) is equal to cII, 
where 11 is the (unsigned) area within the curves given by 
the image of the solution, the set Cq and the set V (hatched 
area in Figure |^. 

The monotonicity of the dissipation is clear from Figure 

As a consequence, there is only one initial condition on 
Cq for which (|^ holds. 

Theorem 1: Under Assumption [T] there exists a unique 
nontrivial hybrid periodic orbit for the hybrid system (|^. 
Tb{x) = Tf{x) at each point x of the hybrid periodic orbit. 



Fig. 4; The hatched areas is proportional to the energy 
dissipated by damping. 

IV. Global asymptotic stability 
The stability of the nontrivial hybrid periodic orbit is a set 
stability problem. Consider the attractor given by 

A = {x&C:Tb{x)=Tj{x),xfQ}. (7) 

Energy considerations similar to those in the previous section 
readily show that A is compact and forward invariant (see 
the proof of Lemma [^. The images of all nontrivial hybrid 
periodic trajectories of Q coincide with A. Convergence and 
stability of the periodic motion follow from the next theorem. 

Theorem 2: Under Assumption [T] the set A in ([7l is an 
asymptotically stable attractor for the hybrid system with 
basin of attraction Ba = M^\{0}. 


























The origin a; = 0 is not in Ba because it is a weak 
equilibrium; solutions to Q starting from the origin can flow 
forever staying at the origin or may jump to —0 or 0 and 
then converge to the hybrid periodic orbit. 

We remark that the stability of the set A does not require 
an explicit characterization of the hybrid periodic orbit. 
We only need to ensure the feasibility of the balance in 
Therefore, by Theorem [T] and Theorem the reset 
feedback law induces a hybrid periodic trajectory for every 
parameters selection that satisfies Assumption [T] Then future 
work comprises performing optimal selections of 0 (and the 
arising periodic motion) via hybrid adaptation. 

The proof of TheoremSis based on a Lyapunov argument. 
Using the definitions in Q, consider the Lyapunov function 
candidate for A given by 

, {u,{x) + n{x)-Tf{x)-\ke^f 

= ( 8 ) 
Ub{x) 

The shape and the level sets of V are illustrated in Figure 
Note that the Lyapunov function V blows up as x approaches 
0 (the boundary of Ba) and as x grows unbounded, as shown 
in Figure for the same parameter selection of Figure 



Lyapunov function V(x) 


Fig. 5; The Lyapunov function V in logarithmic scale. 

The next lemma is a key step for proving Theorem!^ 
Lemma 2: Under Assumption [T] the set ,4 in Q is 
nonempty and compact and the Lyapunov function U in ([^ 


(i) is positive definite with respect to A on Ba n{CU'D), 
namely 

V{x) = 0 if X G A 
V{x) > 0 if a; G n (C U V)\A 
lim V ix) = +00 

|x|->-0+ 

|a:|—>-+oo 

(ii) is constant in the flow directiorj^ 

{VVix)Jix))=0,yxeC (9b) 

^We do not give a formal proof of smoothness of V in our derivation, 
therefore it would be more appropriate to use the directional derivative of 
V in j9b^ . We use here the notation with the gradient to keep the discussion 
simple. 


(iii) provides strict decrease across jump^ 

V (G(x)) - V(x) <0, VxGBAn V\A. (9c) 

Remark 2: Following [1, Corollary 7.32] and [1, Defini¬ 
tion 7.29]), item of Lemmaj^implies that for any indicator 
function w of .4 on Ba there exists class /Coo functions a 
and a such that 


a(tLi(x)) < U(a:) < a{uj{x)), (10) 

which entails standard additional features of A, e.g., ro¬ 
bustness and semiglobal-practical robustness of asymptotic 
stability of A [1, Chap. 6.7]. j 

Based on Lemmathe proof of Theorem]^ is a mere ap¬ 
plication of fundamental Lyapunov results holding for hybrid 
systems described with the formalism in [1]. In particular. 
Lemma establishes that function U in ([^ is a non-strict 
Lyapunov function for compact attractor A, in the sense of 
[11, Equations (6a), (7)]. Indeed, the fact that A is compact 
is sufficient to obtain that ( |9a) i implies ( |T0) | for any indicator 
function ui and suitable class /Coo functions a, a. Moreover, 
(j9^ coincides with [11, Equation (7a)] and ( [9^ implies [11, 
Equation (7b)] for a suitable positive definite function p, once 
again because A is compact. As a consequence, we may 
apply a local generalization of [11, Theorem 2] by noticing 
that all solutions to 0 are semiglobally persistently jumping, 
namely for any arbitrarily large number A, we have that 
solutions restricted to flo'wj^ in C n AB n Ba and jump from 
V n AB n Ba have a uniform reverse dwell time (i.e., a 
maximum time between each pair of consecutive jumps in 
their domain). Such a maximum time is easily computed 
as the minimum flow time from the vertical portion of Cq 
intersected with AB. Note, in particular, that by homogeneity 
of the linear flow equation, such a maximum flow time is 
strictly smaller than the flow time of the first flowing interval 
of the (unique) solution starting from x(0) = [^]- 

V. Simulations 

We simulate the hybrid system 0 with the same param¬ 
eters adopted in Eigure corresponding to m = 1 kg, c = 
0.3 Ns/m, fc = 1 N/m (eigenvalues si 2 = —0.15 ± /0.9887, 
consistent with Assumption [^. Instead of 0 = 0.2 m we 
enforce a reset to the larger value 0 = 0.3. Compared to 
Eigure the hybrid periodic orbit changes as shown in the 
upper part of Eigure]^ The bottom part of Eigure]^ shows the 
values of the Lyapunov function and of the states {xi,X 2 ) 
as functions of ordinary time. 


VI. Conclusions and future research 


Eollowing [6], different hybrid periodic orbits can be gen¬ 
erated by exploiting different reset laws. Within the hybrid 
characterization 0^ one of the switching laws proposed in 
[6] is captured by the flow and jump sets C and T> given in 


Eigure 


7a 


> 2e^ and > 0 is a new control parameter. 


^Note that sgn(a; 2 ) is single-valued for all a: S D n Ba = {x £ 
: xi = 0, X 2 A 0} because it is set-valued only at X 2 = 0. 

^As customary, we denote by B the closed unit ball. 
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Fig. 6: Upper part; phase plot of the trajectories. Bottom 
part: Lyapunov function and hybrid trajectories projected 
on ordinary time for the two initial conditions. Left; initial 
condition (0.1,—0.05). Right; initial condition (0.5,—0.05). 


G{x) = a:i+e^ign(£C 2 ) and f{x) remains ( |^ with the 
parameters of Section The simulation of a number of 
trajectories is reported in Figure [7b] The approach and results 
in Sections|I^and|^can be extended to capture the stability 
properties of this new hybrid system, to show the existence of 
a globally asymptotically stable hybrid periodic orbit. Figure 



(a) Flow and jump sets. 


Phase plane 



(b) Trajectories. 


Fig. 7: = 0.2 and 6 = 0.6. 


[^ shows the potential of the hybrid characterization in (j^. A 
number of reset feedback laws can be modeled by variations 
on the definitions of the sets C and V and of the jump map G. 
The monotonicity of the dissipated energy is crucial for the 
uniqueness of the attractor. Future research will investigate 
the family of reset feedback that guarantees existence and 
uniqueness of an asymptotically stable hybrid periodic orbit. 

It is natural to consider mechanical systems with elastic 
potentials typical of nonlinear springs, so that ( j^ would 
read ±2 = — ^ where the potential U is any posi¬ 

tive definite function. A necessary condition for the existence 
of an attractor with basin K^\{0} is the strict monotonicity 
of the elastic potential, namely the positivity of the function 
tci I—> xi ■ In fact, the lack of strict monotonicity may 

lead to the coexistence of several attractors. The minimal 
requirement for the existence of globally attractive hybrid 


periodic orbits is that the linearization of the flow dynamics 
at the origin has complex conjugate eigenvalues, which is a 
straightforward generalization of Assumption [T] Otherwise, 
trajectories that are sufficiently close to the origin would flow 
towards it remaining in the second/fourth quadrant, without 
triggering any reset. Future work will look for tight sufficient 
conditions for the existence of a global attractor for general 
elastic potentials. 

Finally, the generality of a Lyapunov approach for stability 
analysis calls for extensions of the method to more general 
mechanical systems. The first step is the analysis of n- 
dimensional linear mechanical systems perturbed by resets. 
We also recall that the model in (j^, albeit an abstraction, 
stems from the robotics context. Further detailed analysis 
will investigate its potential for robotic applications, in 
particular the one of legged locomotion. 

Appendix 


A. Proof of Lemma [^ 

The work performed by the nonconservative viscous force 
Fd = cx 2 in moving the point mass from position 1 to 
position 2 causes a change E 2 — Ei in the total mechanical 
energy [10, Page 9]. With the aid of Figurej^ let us exemplify 
it on trajectories like 1 and denote xim the value at which 
a trajectory crosses the line X 2 = 0, so that on the trajectory 
X 2 can be expressed as function of xi for each half-plane 
X 2 > Q, X 2 < 0. Then, by splitting the integral relative 
to the work in two pieces, we get E{x{t 2 )) — E{x{ti)) = 
ft'^-(cx2}dxi + J°^^-(cx2)dxi = -cU. 


B. Proof of Lemma [^ 

The proof of this fact relies on Lemma From (j^, 

EbG) Ef{x) 

Ub{x) + Tb{x) - Tf{x) - 
Ub{x) 

(cn(a:) - Uf 
U^) ’ 

where Eb and Ef are the total energies of the system right 
after and right before a jump, respectively; t/ is the potential 
energy at [®] S Co; n(x) is the area spanned by the solution 
passing through x during a flow from Cq to V (where Ef 
is evaluated); and c > 0 is the damping coefficient in Q. 
Figure 1^ provides two examples for n(a:). 

We are now ready to prove Lemma [^ First notice that, 
due to the uniqueness of flowing solutions, the function 
X n(a:) is necessarily strictly increasing as x moves 
farther from the origin (or any compact set). Denote by 
cIIo = cll ([ Q ]) the dissipated energy when starting from 
the corner of set Cq in 0- Moreover, denote by cII* = 
= U the dissipated energy associated to the hybrid 
periodic orbit. Note that 11* > Hq necessarily, because cIIo 
cannot be larger than the total energy U = cII* at the 















































































beginning of the corresponding solution starting from the 
corner of Cq. Then, 

A = {xGC: Il{x) =n\x^ 0}, (12) 


which proves that it is non-empty and compact. We prove 
now the three items of the Lemma. 

Item (j^. Since V(x) in (111 is non-negative and zero if 
and only if n(a;) = If*, from expression ( |T2| ) we obtain 
V(x) = 0 if and only \f x A and positive otherwise. 
Moreover, as x approaches zero, we have that Ub{x) tends 
to zero, which implies V{x) —> oo. Since Ub{x) < U for all 
X due to the structure of Cq, lim| 2 ;|_>,oo n(a;) = -foo implies 
lim|x|->oo ”1^(2;) = +00. 

Item (j^. This item follows in a straightforward way from 
the fact that V(x) remains constant by construction during 
flow. 

Item ([ilill. First of all notice that x = [%l] & 'D implies 
esign(a;2) for all X S V n Ba- 
We split the proof in three cases. We only consider jumps 
from points in the negative part of V (namely x -2 < 0) 
because of the central symmetry of the phase portrait. We 
also use the simplified notation 11+ to denote n(x+). Similar 
simplifications will be used for other quantities. 

Case 1: If > If* > Hq. First of all, by uniqueness of 
solutions n+ > If*, otherwise the flow would intersect the 
hybrid periodic orbit. Consider the left part of Figure and 
note that If > If* implies 11+ < If. Indeed, exploiting = 
Ub = U = cII* and Tb + Ub = Tf + cII and T^+ = Tf, we 
get Tb = Tf + cIl — cII* > Tf = Tj+. If > 11+ follows by 
monotonicity. Finally, the result is proven from 


xi = 0 and that G(x) = 


0 = cH* - t/ < cn+ - L < cH - (7. 


Case 2: If* > If > IIo. First of all, 11+ < If* again from 
uniqueness of solutions. Consider the right-part of Figure 
and note that IIo < If < If* implies 11+ > If. In fact, 
following the argument of Case 1, we have that = Ub = 
U = cII* and Tb + Ub = Tf + cll and T^ = Tf thus 
Tb = Tf + cli — cII* <Tf = Tj+. The result is proven from 

0 = cn* - U > dl+ - tf > cll - U. 




Fig. 8: Left; If > If* > IIo. Right; 11* > 11 > Hq. Pink 
areas are 11+ mirrored about the origin. 


Case 5.' 0 < n < IIq. Consider Figure]^ In this case we 
have Ub < U because Ub is evaluated on the horizontal part 
of Cq. Then 

V+-V = (cn+ - (7)7(7 - (cH - Uf/Ub 
< (cn+ - uf/u - (cn - uf/u. 


Now, observe that cn+ < U = cII* because otherwise the 
forward solution from x+ would intersect the flowing portion 
of the hybrid periodic orbit (thus contradicting uniqueness). 
Then, using Ub > 0, we get L+ — 1+ < 0 from 

(cn+ - (7)2 - (cn - (7)2 

= c(n+ - n)(cn+ -u + cii -_l) < o, 

>0 <0 <0 

where we used (i) n+ > H from n+ being evaluated from 
the vertical part of Cq and (ii) cH < (7 from cll <Ub. 



Fig. 9; Case 3; 0 < n < Hq. 
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